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ABSTRACT 



Aims. This paper is an H band interferometric examination of Arcturus, a star frequently used as a spatial and spectral calibrator. 
Methods. Using the IOTA 3 telescope interferometer, we performed spectro-interferometric observation (R « 35) of Arcturus. 
Atmospheric models and prescriptions were fitted to the data to derive the brightness distribution of the photosphere. Image recon- 
struction was also obtained using two software algorithms: Wisard and Mira. 
Results. An achromatic power law proved to be a good model of the brightness distribution, with a limb darkening compatible with 
the one derived from atmospheric model simulations using our marcs model. A Rosseland diameter of 21.05 ± 0.21 was derived, 
corresponding to an effective temperature of T^ = 4295 ± 26 K. No companion was detected from the closure phases, with an upper 
limit on the brightness ratio of 8 x 10~ 4 at 1AU. Dynamic range at such distance from the photosphere was established at 1.5 x 10~ 4 
(lcrrms). An upper limit of 1.7 x 10~ 3 was also derived for the level of brightness asymmetries present on the photosphere. 
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1. Introduction three different public ations; iDyck et al.lll996tlPerrin et alJ[i998t 

IVerhoelst et al.ll2005l) . The initial goal of a new observation run 

Arcturus' curse is to be too popular. Many new generations was t0 i everage me more extended capability of IOTA to show 

of instruments - including interferometers - observe it as a the abm of the inter fe rom eter to do a reliable image of a com- 

test object. The reasons for that: this red giant is bright, large, mon i v observed obiect 
and spectrally well-defined. The curse is: since the instruments 

are new, the star usually takes face on the bizarre systematic Indeed, image reconstruction is difficult. Even though it is 
errors of challenging observations. We can cite, among oth- routinely performed by the current generation of radio interfer- 
es inexact diameter and temperature measurements (prompt- ometers, this technique - herein called "regularized imaging" - 
ing lGriffin & Lvnas-Grav ( 1999) to write an article entitled "The remains marginal in optical interferometry. This is simply due 
Effective Temperature of Arcturus"), or false dupli city observa- to a lack of s P atial frequency coverage. Optical interferometers 
tions ("Arcturus as a Double Star" by Grifnn|[l998|). ' are usually more difficult to build, and the complexity quickly 
This paper does not pretend to be an exception - interfer- increases with the number of telescopes. Therefore, since the 
ometry is still a challenging technique. The main difference amount of information accessible in the Fourier plane is sparse, 
is in the interferometer used: at the time of our observations, our ability to reconstruct a reliable image of a complex object is 
IOTA was more wise than new (we shall note that Arcturus nmitea. 
has already been observed several times by IOTA and led to A more common data analysis is to suppose the object 

to be conform to a model - or prescription. Originally, it 

Send offprint requests to: S. Lacour consisted in fitting visibility curves of uniform disks (e.g. 

* Postdoctoral Fellow of the Fund for Scientific Research, Flanders Michelson & Peas ell92Uldi Benedetto & Fov|1 986). With time, 
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Table 1. Calibrators 



Calibrator 


Spectral Type 


UD diameter 


1 HD 120477 

2 HD 125560 

3 HD 129972 


K5.5III 

K3III 

G8.5 III 


4.460 ± 0.050 
1.910 ±0.021 
1.540 ±0.020 


Table 2. Arcturus observing log 


Date 
(UT) 


Interferometer 
Configuration" 


Calibrator 
(Table [Q 



2006 May 1 1 


A15-B05-C10 


1,2,3 


2006 May 12 


A15-B05-C00 


1,2,3 


2006 May 13 


A15-B15-C00 


1,3 


2006 May 14 


A30-B15-C00 


2 


2006 May 16 


A35-B15-C25 


2,3 



Configuration refers to the location in meters of telescopes A, B, C 
on the NE, SE and NE arms respectively. 



-20 
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West direction [MA] 

Fig. 1. u-v coverage. Maximum projected baseline length is 37.7 
meters. The lack of high frequency information in the East-West 
direction is due to the geometry of IOTA. 



it inc luded more complicated models, e.g. limb-darkened disks 
(e.g. Ouirrenbach et al. 1996, Mklll), disks with spots (e.g. 
Young et al. 2000, COAST), disks with a molecular envelope 
(e.g. iPerrin et alj |2004, IOTA), etc. Both methods, regularized 
imaging and model fitting, complement each other. The role of 
regularized imaging is commonly to guide the choice and com- 
plexity of a model. The role of model fitting is to obtain the high- 
est precision results on the parameters of the model. The pitfall 
may be when the model does not best suit the object, hence the 
need for quality regularization imaging. 

Recently the maturation of interferometric facilities (e.g. 
IOTA, CHARA, VLTI) reached the point where u-v cover- 
age (both in amplitude and phase) a llows regularized imaging 
(Monnier et alj|2007at lLacour||2007l) . Here we present the data 
on Arcturus, used as a test star for optical interferometry re- 
construction softwares. Several astrophysical issues also justify 
this investigation: 1) what is the limb dark ening? Is it com- 
patible with red giant atmosphere modeling dDavis e t al. 2000; 
Claret 2000)? 2) are the disputed previou s detections of a com- 
panion compatible with our observations ( Perrvman et all 1997; 
IVerhoelst et al 2005; Brown 20071)? 

The outline of this paper is as follows. Section [2] gives an 
overview of the IOTA interferometer, describes the data reduc- 
tion process and shortly present the dataset. Section[3]compares 
our data with atmosphere models, using either limb-darkening 
prescriptions or a more evolved atmospheric simulation (the 
Marcs model). Section [4] investigates on a possible deviation 
from point symmetry. Results of image reconstruction are pre- 
sented in Section[5j and Section[6]conclude. 



2. Observation and data reduction 

2.1. Description of IOTA observations 

The interferometric data presented herein were obtained us- 
ing the IOTA (Infrared-Optical Telescope Array) interferometer 



dTraub et alJl2003l) . a long baseline interferometer which oper- 
ates at near-infrared wavelengths. It consists of three 0.45 meter 
telescopes movable among 17 stations along two orthogonal lin- 
ear arms. IOTA synthesizes a total aperture size of 35x15 m, cor- 
responding to an angular resolution of « 10 x 23 milliarcseconds 
at 1 .65 jum. Visibility and closure phase measuremen ts were ob- 
tained using the integrated optics combiner IONIC dBerger et al.l 
2003); light from the three telescopes is focused into single- 
mode fibers and injected into the planar integrated optics (IO) 
device. Six IO couplers allow recombinations between each 
pair of telescopes . Fringe detection i s done using a Rockwell 
PICNIC detector (iPedretti et al.ll2004 . The interference fringes 
are temporally-modulated on the detector by scanning piezo mir- 
rors placed in two of the three beams of the interferometer. 

Observations were undertaken in the H band (1.5 fj.m< A < 
1.8 /im) divided into 7 spectral channels. The science target ob- 
servations are interleaved with identical observations of unre- 
solved or partially resolved stars, used to calibrate the inter- 
ferometer's instrumental response and effects of atmospheric 
seeing on the visibility amplitudes. The calibra t or sou rces 
were chosen i n two different catalogs: iBorde et all d2002l) and 
Mer and et alj d2006l) : using criteria on the separation (< 10 de- 
grees) and magnitude. The calibrators are listed in Table Q] 

Arcturus was observed in May 2006 during 5 nights and us- 
ing 5 different configurations of the interferometer. Full observa- 
tion information can be found in Table[2] including dates of ob- 
servation, interferometer configurations and calibrators. Fig. Q] 
shows the u-v coverage achieved during this observation run. 
The geometry of the IOTA interferometer and the position of the 
star on the sky constrained the extent of frequency coverage. We 
covered a frequency range equivalent to the one of an elliptical 
telescope of aperture 38x15 meters, with a 20 degree inclination 
East of North. 

2.2. Data reduction 

Reduction of the IONIC visibility data was carried out using 
custom sof tware similar in its main princ iples to the one de- 
scribed by ICoude Du Foresto et al.l dl997). We measured the 
power spectrum of each interferogram (proportional to the tar- 
get squared visibility, V 2 ), after correcting for intensity fluctua- 
tions and subtracting bias terms fro m read noise , residual inten- 
sity fluctuations, and photon noise (Perrin 2003). Next, the data 
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Fig. 2. Single scan obtained on calibrator HD 120477. It shows the flux (arbitrary units) as a function of the optical path difference 
(OPD) for the three baselines, and the seven spectral channels. Eye inspection allows to observe the decreasing frequency - and 
therefore increasing wavelength - from top to bottom. The bottom fringes correspond to the sum of the spectral channels, showing 
a radical change in the coherence length. Each datapoint is composed of 200 scans. 



pipeline applies a correction for the variabl e flux ratio s for each 
baseline by using a flux transfer matrix (Monnier 2001). Finally, 
raw squared visibilities are calibrated using the raw visibilities 
obtained by the same means on the calibrator stars. Calibration 
accuracy had been studied by extensive observation of the b inary 
star A Vir. For bright stars (H mag < 5), IZhao et all (|2007) have 
validated 2 % calibration error for V 2 ; corresponding to a 1% 
error in visibility. We therefore systematically added a 2% cali- 
bration error on all the squared visibilities present in this paper. 
In order to measure the closure phase (CP), a fringe track- 
ing algorit hm was applied in real-time while recording interfer- 
ograms (Pedre tti et alj|20 05), ensuring that interference occurs 
simultaneously for all baselines. We required that interferograms 
are detected for at least two of the three baselines in order to as- 
sure a good closure phase measurement. This technique, called 
"baseline bootstrapping" allowed precise visibility and closure 
phase measurements for a third baseline with very small coher- 
ence fringes. We followed the method of Bald win et alj 0996) 
for calculating the complex triple amplitude and deriving the 
closure phase. Pair-wise combiners (such as IONIC) can have a 
large instrumental offset for the closure phase which requires to 
be calibrated by the closure phase of the calibrator stars. We no- 
ticed very stable closure phase measurements during the nights 
with drifts of less than a degree. 



2.3. Wavelength calibration 

Spectral information was obtained by the means of a prism 
place d between the int egrated optics and the PICNIC cam- 
era (|Ragland et al. 2003). The temporally-modulated fringes are 
therefore spatially dispersed on the detector. To ensure well- 
defined spectral edges, we also inserted a broad band H filter 
in the optical path. Its bandpass is spatially equivalent to seven 
pixels on the camera. 

Wavelength calibration of the spectral channels is however a 
difficult and critical step. This is especially important since the 
prism was removed and re-inserted (with a slightly different po- 
sition) between the night of the 13th and the 14th. Fortunately, 
the spectral wavelength is coded in the data (see Fig. |2). The 
fringe frequency (in pixels 1 ) is directly proportional to the 
wavenumber. The factor of proportionality is constant since the 




E 



X [urn] 

Fig. 3. Relative photometry between the different spectral chan- 
nels as a function of the wavelength (arbitrary vertical units). 
The wavelength was determined by measuring the frequency of 
the fringes as shown in Fig. [2] 



modulation of the optical path is done by moving the piezo mir- 
rors a certain distance (step-like) between each pixels reads, even 
though the steps are smoothed out by mirror/mount inertia. The 
relative wavelength between each channel and each night was es- 
tablished this way with a precision better than 0.1%. This level 
of precision was achievable thanks to small differential piston 
variations due to good seeing conditions and fast reading mode. 
The relative wavelength between each baseline was also studied. 
To do so, we compared the fringe frequency observed at a given 
spectral channel between the three baselines. The frequency of 
the third baseline is equal to the sum of the frequency of the 
two first, within 0.2% error bars. This means that the different 
baselines are at equal wavelength at a 0.2% level. 

However, absolute calibration requires to know the exact an- 
gle of the incoming beam on the piezo mirror. March 2007 nar- 
row band observations were used, and allowed to establish the 
speed of the optical path modulation at 0. 1 88 ± 0.002 yum/sample 
for the first delay line, and 0.195 ± 0.002 //m/s ample for the sec- 
ond delay line. Optical path modulation was measured on the 
third baseline at 0.383 + 0.004 /im/sample. A time sample cor- 
respond to the integration time between two reads. The » 1% 
error bars are mainly due to uncertain changes in the angle of 
reflection which occurred between March and May 2006. 

Figure [3] summarizes the wavelength calibration results by 
plotting the average integrated flux on each spectral channel as a 
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Fig. 4. Overview of the dataset: Visibility square measurements 
as a function of the baseline length. The wavelength is color 
coded with the same palette used in Fig. [I] The solid curve corre- 
spond to the visibility curve of a uniform stellar disk of angular 
diameter 20.30 mas (not accounting for bandwidth smearing). 
The bottom panel presents the residual of that fit, showing the 
clear inconsistency of the second lobe. 



















CO 

CD 
CD 

O) 
CD 

2, 


100 












" 


CD 
CO 











o 
O 


y 




CC 
Q. 


C10 


o 
O 


o 
O 




LO 

CJ- 

O 


CD 

3 




LO 

m 


LO 

m 


LO 

m 


LO 

m 




LO 


CO 

O 

o 


-100 


A15- 


A15 


LO 

< 


A30- 




LO 

co_ 

< 



100 



200 



300 



CD 

CO 

CO 

-C 
Q. 
CD 

CO 

o 



en -10 r 



CD 




100 



200 



300 



Observation number 



Fig. 5. Overview of the dataset: closure phase measurement as a 
function of observation datafile number. The bottom panel shows 
the residual of fitting a simple limb darkening disk to the data 
(power law with parameters as stated in Table [3}. Color legend 
is similar to the one used in Fig.[T]and[4] Because such a model 
corresponds to a symmetrical brightness distribution, the closure 
phases are either or 1 80 degrees. The reduced^ 2 of the closure 
phase alone is 1.06. 



function of wavelength. The fairly large error bars in wavelength 
are mainly due to the uncertainty on the angle of reflection, and 
correspond to a possible global shift in absolute calibration. In 
contrast, relative wavelength is precisely established, and shows 
that a significant displacement of the prism occurred between the 
night of the 13th (squares), and the night of the 14th (cross). 

2.4. IOTA field of view 

The high resolution of IOTA has a counterpart: the field of view 
is limited. The first limitation is due to the field of view of each 
individual telescopes, delimited by the cone of acceptance of the 
fibers on the sky. Such value is difficult to estimate, since it de- 
pends on the interferometer as well as the atmospheric seeing. 
A first order estimation is to neglect the atmospheric turbulence 
and to consider the fiber core to be filling the diffraction pat- 
tern of the telescopes. In this assumption, the field of view of the 
telescopes reads: 



FOV, 



telescopes — 



D 



(1) 



where D is the diameter of an individual telescope. 

A second limitation is the field of view of the interferometer. 
It is delimited by the maximum distance between two objects 
which fringes overlap on the detector. To be rigorous, one should 
take into account parameters like the mode of recombinaison, the 
stroke of the piezo (in the case of IONIC), and even the spectral 
energy distribution of the target. However, to establish a simple 
relation, we will only take into account the spectral bandwidth 



of a spectral channel (AA), as well as the distance between two 
telescopes (B): 



FOVi 



interferometer — 



AAB 



(2) 



Note that the interferometric field of view is baseline dependent. 
It will be larger for shorter baselines, and smaller for longer 
baselines. Moreover, this field limitation is valid only in the 
direction along the baseline. Perpendicular to the baseline, the 
bandpass does not cause any field limitation. All in all, it is diffi- 
cult to establish the field of view of an interferometer as a whole. 
A conservative way to do so is to consider the maximum baseline 
length for a given direction. 

In the North-East/South- West direction, using the spectral 
dispersion mode of IOTA (D - 45 cm, A/l = 40 nm and 
B — 35 m), the field of view is not limited by the tele- 
scope (.FOVteiescopes = 750 mas), but by the bandwidth. The 
field of view is 350 mas at 1.55^/m, and 480 mas at 1.80//m. 
In the North-West/South-East direction, the shorter baselines 
(B - 15 m) allow a larger interferometric field of view, hence 
a 750 mas field limitation due to the telescopes size. 

We will consider in the following of this paper a 
400x750 mas field of view for IOTAB 



IOTA's field of view decrease when using large band filters 
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Table 3. Diameter and limb darkening measurements 



Law 


parameters 


Reduced^ 2 


Uniform 


<9 UD = 20.304 ±0.011 mas 


31 


Power 


<9 LD = 20.900 ± 0.007 mas 
a = 0.258 ± 0.003 


1.962 


Quadratic 


9 LD = 20.922 ±0.036 mas 
a = 0.186 ±0.021 
b = 0.298 ± 0.053 


1.959 


Quadratic 
Claret (2000) 


9 LD = 20.931 ±0.004 mas 
a = 0.0291 
b = 0.5107 


2.956 



Non-linear 
Claret (2000) 



6 LD = 20.863 i 
ai = 0.8175 
a 2 = 0.0827 
a 3 = -0.4116 
a 4 = 0.1864 



0.004 mas 



2.013 



marcs model ft 



Ross 



21.05 ±0.01 mas 



2.080 



Note - Errors bars are pure calculations based on the second derivate 
of the^ 2 . They are not valid when assuming an unrealistic model of 
the brightness distribution (for example a uniform disk). Moreover, 
diameter errors do not include the 1% uncertainty due to an eventual 
wavelength miscalibration (see section l2~3l l. 



2.5. The dataset 

The dataset consist of 924 visibility measurements and 308 clo- 
sure phases. The V 2 are plotted as a function of the baseline 
length in the upper panel of Fig. [4] The CP are plotted in the 
upper panel of Fig. [5] The frequency plane coverage was previ- 
ously presented in Fig.Q] The solid curve corresponds to the best 
fit of a uniform disk. The residuals are plotted on the lower pan- 
els. It is interesting to note that fringes have been observed with 
a contrast below 1 %. Such a low contrast exists thanks to the dis- 
persive mode, which allows a deep first null. If the full H band 
was observed, the effect of bandwidth smeari ng would have lim- 
ited th e depth of the null to several percents dPerrin & Ridgwavl 
2005). Probing into the null was possible thanks to bootstrap- 
ping, two baselines of sufficient contrast being enough to track 
the fringes on all the baselines. 

A few things are striking: first, the object is relatively achro- 
matic. This can be seen on the residuals of the V 2 . Secondly, the 
second lobe of the data is not well fitted by a uniform disk. This 
is due to the presence of limb darkening. Thirdly, the closure 
phases are close to zero or n. It means the object is likely to be 
point symmetric. 

3. Comparison with atmosphere models / 
prescriptions 

3.1. Fitting limb-darkening prescriptions 

Since limb darkening is apparent, a logical first step is to 
fit a model for the brightness distribution of the photosphere. 
Numerous types of limb-darkening (LD) prescriptions exist in 
the literature. We us ed two of them, which we supposed achro- 
matic. A power law dHestrofferll 19971) : 



and a quadratic law dManduca et al]fl977l) : 



(3) 



(4) 



where // = -yl - (2r/#LD) 2 X r being the angular distance from 
the star center, and #ld the angular diameter of the photosphere. 



In terms of complex visibilities, the power law limb darkening 
prescription yields: 



^r(Q</2 + 



T(a/2 + 2) 



r(a/2 + k + 2)r(k+l) 



-(nv r LD y 



(5) 



where v r is the radial spatial frequency and Y the Euler function 
( T(k + 1) = k\). On the other hand, the quadratic law yields: 



n uMO . a + 2bh l2 {Q J 2 (0 

(1 - a - b) — : 1 — -tt; 2b- 



V(v r ) = 



1/2 -a/6 -ft/12 



(6) 



where £ = 7tv,0ld, J\ and J 2 are the first and second-order Bessel 
function respectively, and: 



'^O-^f^-cosu.j. 



(7) 



Results for the fits are presented in Table[3] Using Eq. (0), we 
obtained a^ 2 of 2413, for 1230 degrees of freedom. The reduced 
X 2 (x 1 over the number of degrees of freedom) does not improve 
significantly when using a two-parameter prescription for limb 
darkening, prompting us to consider the power law model as a 
sufficient approximation. 

There are two main explanations for the reduced x 1 differ- 
ent from one: (i) an underestimation of the error bars, and (ii) 
an inexact prescription of the brightness distribution by Eqs. ([3]l 
and (|4]i. Making the distinction between these two is difficult. 
On the one hand, photometric variations of the star are observed 
at the order of one percent ( Retter et al.l 2003). an indication that 
the brightness distribution should not be as smooth and symmet- 
ric as our prescriptions are. On the other hand, no deviation from 
point symmetry is observed in the closure phases (section |4| 
whose reduced^- 2 , taken independently, is 1.015. The depar- 
ture from simple LD models could therefore only be explained 
by a missing point-symmetric component. The residuals are dis- 
cussed more throughly in the last paragraph of Section [3~2l 

Whatever the cause, we decided to be as conservative as pos- 
sible by scaling the errors to a x 2 of one. The error bars stated 
in Table [3] are obtained by this mean. To decrease the x 2 , we 
explored - and discarded - two other alternatives. The first one 
was to increase the error due to calibration (higher than the 2% 
justified in section l2~2l . However, this dramatically increased the 
errors bar on the lowest frequencies, which is not desired since 
they are already well fitted by our prescriptions. The second ap- 
proach was to add an additive error due to a potentially imper- 
fect subtraction of the power spectrum bias. Such an error was 
included at a 2% level for faint objects by Monnier et al. (2006). 
However, this dramatically increased the errors bars on the low- 
est visibilities, which brings an unnecessary bias on data whose 
first zero is already well fitted. In conclusion, a global scaling of 
the error bars was seen as the best alternative. 



3.2. Fitting the Marcs model 

The M arcs atmospheric model was presented in lVerhoelst et al.l 
(2005). The models were originally constructed and fine-tuned 
for the calibration of the ISO-S WS (Infrared Space Observatory 
Short Wavelength Spectrometer) and checked against FTS spec- 
tra (lDecin||2000h . For the present study, we searched the full 
Arcturus FTS spectral atlas (Hinkleetal. 1995) in the H band 
for peculiar spectral features. Lines are sparse and well spaced. 
They belong mainly to CN, OH and some atomic transitions. The 
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Fig. 6. The synthetic H-band spectrum of the marcs model (solid 
line: in opacity sampling resolution, dotted line: convolved to the 
instrumental spectral resolution - shifted up by five thousands 
Jensky) and the central wavelengths of the spectral channels of 
IONIC. The peak in the spectrum corresponds to the H opacity 
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Fig. 7. Wavelength dependent visibility curves derived from the 
Marcs atmospheric model. The color coding is similar to the 
one used in Fig. [6] IONIC data are superimposed on the curves, 
and residuals are plotted in the lower panel. The closure phase 
residuals are identical to those of Fig. [5] 



IONIC data are therefore ideal to study the H continuum, which 
has its minimum (the transition between bound-free and free- 
free regimes) within the bandpass sampled by our data. The only 
free parameter to match model to observations is the angular di- 
ameter corresponding to the outermost point in our model inten- 
sity profiles (tr oss = 10~ 7 ). Several models with stellar parame- 
ters around those determined by Decin et al. (2003) were used, 
but they bring no significant improvement in^- 2 compared to the 
spectroscopically preferred model (T e ff = 4320 K, logg - 1.5, 
[Fe/H] = -0.5 and v turb = 2 km s" 1 ). 

The synthetic H-band spectrum calculated from our model 
and the comparison of our dataset with this model are shown 
in Fig. [6] and Fig. [7] We find the best agreement for a diameter 
of 21.465 ± 0.008 mas, which corresponds to a tr oss = 1 diam- 



eter of 21.05 mas. This diameter is slightly larger than the one 
found in Sect. 13. II the star appears a little smaller at wavelengths 
of minimal photospheric opacity than at the Rosseland-averaged 
opacity. 

With ax 2 of 2, this fit is almost as good as it was when fitting 
a free-parameter limb darkening prescription. This is an overall 
confirmation of the validity of the Marcs modeling of the limb 
darkening. Analysis of the visibility residuals plotted in Fig. [7] 
indicates some possible shortcomings of the model (supposing 
error bars are not underestimated, see Section [3~TT i. Indeed, the 
high^ 2 can be accounted for by two biases: a chromatic bias at 
low frequency (« 9 MA), another achromatic at high frequencies 
(around the second nul). Accounting for these biases could be 
done by (i) introducing a circumstellar emission of H 2O at a level 
of ha lf a percent (water detection was reported by Rvd e et al.l 
2002, although normal hydrostatic model does not predict any in 
the photosphere), and (ii) slightly modifying the limb darkening 
distribution. However, such possibilities are at the limit of what 
we think is reasonable to derive from our data, and no further 
modeling was done to avoid over-interpretations. 

3.3. On the angular diameter of Arcturus 

Numerous angular diameter measurements can be found in 
the literature. Previous interferometric observations either use 
uniform-disk fitting and apply limb-darkened corrections, or fit 
disks whose limb-darkenin g is fixed by atmospheri c models. 
At a wavelength of 2.2yum, di Benedetto & Foy ( 1986) observed 
Arcturus with the I2T interferometer and published a diame- 
ter of #ud = 20.36 ± 0.20 mas as well as a limb-darkened 
value #ld = 20.95 ± 0.20 mas. Previous measurements using 
the IOTA interferometer exist too, and yi elded in the K band 
6» LD = 19.5 ± 1.0 mas (iDycket all 119961 usi ng bulk optics), 
6 LD = 20.91 ± 0.08 mas "(Perrinet ailll998l u sing FLUOR) 
and #r oss — 21.18 + 0.21 mas (Verhoelst et al. 2005, also us- 
ing FLUOR). In the visible. iMozurkewich et al.l (12003) observed 
Arcturus using the Marklll interferometer (450-800 nm), and af- 
ter correction for a substantial limb darkening effect, published 
6 hD = 21.373 ±0.247 mas. 

From our dataset, and taking into account wavelength cali- 
bration uncertainties, we derived Old = 20.91 + 0.21 mas and 
#Ross = 21.05 + 0.21 mas. These results are in agreement wit h 
I2T, Marklll and IOTA observations (1.5cr in lDvcket all 19961) . 
It does not yield an increase in terms of precision, but our mea- 
surements are indeed interesting since, unlike the others, they 
did not require a pre-defined value to a ccount for limb dark- 
ening. Using the Rosseland diameter and Griffi n & Lvnas-Gravl 
(1 19991) estimation of the integrated flux (F = (4.98 ± 0.02) x 
10 _5 ergcm _1 s -1 ), we can update their calculation of Arcturus' 
effective temperature to r e ff = 4295 + 26 K. 



3.4. On the limb darkening 

An important prospect of this work was to compare our limb 
darkening measurements with existing atmospheric models. A 
first test was to derive parameters of the limb darkening, and 
compare them with published values in the literature. We were 
greatly surprised to see a strong inadequacy bet ween our mea- 
surements and the quadratic parameters given by IClaretl ([2000) 
(see Tableland Fig. [8}. However, they also published the values 
for a more complex 4-parameter non-linear law: 



w 

1(1) 



= l-ai(l- ^D- fl2 (l-/z)-a 3 (l V /2 )-a 4 (l-^).(8) 
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Fig. 8. Intensity profiles of our marcs model as colored lines. 
For comparison, the best parametric fit of a power law is repre- 
sented with diamonds, the Kurucz model with triangle, and the 
quadratic LD curve IClaretl (|2000) as a dotted line. This last fit 
differs significantly form the others, revealing a problem in the 
method used for limb darkening fitting used by Claret. 
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Fig. 9. Limb darkening as a function of the wavelength. Upper 
panel: the observed LD coefficients a (with error bars) are com- 
pared with the LD coefficients derived by fitting a I(jj) - p a 
profile to the Marcs intensity profiles of our dedicated Arcturus 
model (T eff = 4320 K, logg = 1.5, [Fe/H] = -0.5 and v turb = 
2kms~'). Lower panel: influence of temperature and gravity on 
the limb darkening. 



IClaretl J2000h claims that this four-parameter non-linear law 
should give a more reliable estimation of the limb darken- 
ing. Using his published parameters (assuming T e fr = 4250 K, 
logg = 1.5, [Fe/H] = -0.5 and v tur b = 2kms~ 1 ) we were able to 
indeed confirm a correct fit. 

However, it is not the quadratic law which is intrinsically less 
able to match the limb darkening: when leaving the parameters 
free to adjust, the x 1 of a quadratic law is able to match the x 2 
of the non-linear law . Therefore, th e problem with the quadratic 
values published by Claret (2000) should lie in the method to 
derive the parameters. To confirm this, we used the Atlas model 
(Kurucz 1979) - the one used by Claret - and we were able to 
obtain a good fit for the limb darkening (see Fig. [8}. An expla- 
nation could be found in a paper written by iHevrovskvl (12007). 



in which the author states that conventional stellar limb fitting 
methods (like the one used by Claret) are biased. 

But the most striking results from Table [3]is the consistency 
in fitting quality achieved when using different limb darkening 
laws (except when using Claret's quadratic value). We noted 
that both Marcs and Atlas models give similar fits, show- 
ing an equivalent capacity to correctly model the atmosphere 
of Arcturus. The reduced^- 2 values are not exactly 1, but are 
close to the ones obtained when fitting LD laws with freely vari- 
able parameters. This is a good validation of both atmospheric 
modeling softwares. Secondly, we do not note any difference in 
the fitting quality between a power and a quadratic limb dark- 
ening law. Furthermore, the likelihood does not increase when 
using a 4-coefficient non-linear law (reduced^ 2 of 1.97 - we do 
not present the results in Table [3] since none of the coefficients 
are properly constrained by our dataset). This is because we do 
not have the necessary angular resolution to distinguish the sev- 
eral limb darkening laws used here. To our dataset, all of them 
are equally good. Therefore, for angular resolution no greater 
than ours, we recommend using the power law instead of the 
two other limb darkening laws tested in this work, since it would 
use less free parameters while still being able to correctly model 
theLD. 

Finally, we investigated the spectral dependence of the limb 
darkening. To do so, we fitted a limb-darkening power law with 
a wavelength-dependent a value to the observations. The fit was 
done using an achromatic diameter as, in principle, there is no 
such thing as a different diameter at different wavelengths: a dif- 
ferent intensity profile, or in the case of a well-behaved star just 
a different LD, mimics a different diameter at different wave- 
lengths. In the absence of extended molecular layers and other 
similar large deviations from a normal photospheric IP, this ef- 
fect is mostly accounted for by the LD parameter. Similarly, we 
derived theoretical a values from our preferred marcs model. 
The result is summarized in Fig. [9] The general agreement is 
quite good. The almost linear slope of the limb-darkening is in 
fact a complex combination of opacity due to the H contin- 
uum and molecular absorptions features. A minor discrepancy 
is an overestimation of the LD at the blue end of the bandpass. 
We searched a grid of models for possible improvement, but no 
significantly better fit could be attained with reasonable stellar 
parameters. 



4. On the point symmetry of Arcturus 

4.1. Fitting the closure phases 

Closure phases (CP) are extremely sensitive to deviations from 
point symmetric brightness distributions. For example, a binary 
of contrast ratio 1:100 could induce closure phases of several 
tens of degrees at low visibilities. The mean of our first 147 
CP measurements (first two days of observation) is 0.067 de- 
gree, with an average root mean square of 0.34 degree. Such 
high quality data is therefore excellent for probing a companion. 
When fitting a power law limb-darkened disk to the data (both 
V 2 and CP; see section [3~TT >. the^- 2 on the CP was 327 over 308 
closure phases — corresponding to a reduced^ 2 of 1.07. Fig. 
shows the CP as well as the residual of the fit. The fit is, in our 
opinion, satisfactory. 

An upper limit for the brightness ratio of a possible com- 
panion can be obtained. We modified the visibility function pre- 
sented in Eq. (0 to account for the presence of a point-like off- 
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Fig. 10. Non-detection of a companion to Arcturus. Left panel: closure phase x 1 map. Right panel: 1 cr upper limits of the brightness 
ratio of a companion. The circles represent the limit of the photosphere. The field of view of IOTA is 400x750 mas (see discussion 
in Section l2~4l >. 



centered source: 
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-n 2 6 2 (u 2 + v 2 )^ 



+Kexp(2i7r(Xu-Yv)) 



(9) 



K is the brightness ratio of the companion, X and Y its position 
and u and v the spatial frequencies (arcsec 1 ). The star parame- 
ters (a and #ld) are fixed to the value presented in Table [3] For 
each position of the companion Eq. (O is computed, CP are de- 
rived and K is adjusted to minimize the^f 2 on the closure phases. 
The minimum^ 2 are plotted as a function of X and Y on the left 
panel of Fig.[l0] The general minimum^ 2 for a companion situ- 
ated within the field of view of IOTA but further away than 1 AU 
of the star (400 mas> Vx 2 + Y 2 > 89 mas) is 299, with a bright- 
ness ratio K — (4 + 4) x 10~ 4 . This is not significant enough to be 
considered a detection. The values (K + cr(K)) can nevertheless 
be used to derive upper limits for the brightness ratio of a possi- 
ble binary system. It is plotted in the right panel of Fig. \W\ The 
average dynamic range at 1 AU of the star is 1 .5 x 10~ 4 . 

Closer to the photosphere, the x 2 can decrease substantially. 
The minimum is 257, for X — 10 and 7=11 mas. It still does 
not mean we have detected anything, since this value is below 
the number of degrees of freedom. However, the fit can be used 
to put upper limits on the brightness of a possible hotspot on the 
photosphere. The maximum value for K + o~(K) on the photo- 
sphere is 1.7 x 10~ 3 . Note that the signature of an hotspot on the 
CP gets smaller when it gets closer to the photocentre. Therefore, 
we cannot exclude the presence of a bright hotspot coinciden- 
tally situated in the middle of the photosphere. 



keep using Arcturus as a calibrator. Our results put an upper limit 
on the brightness ratio of a possible companion to 8 X 10~ 4 in the 
H band. 

To make our results compatible with a bi nar y system as pro- 

S iosed by Hipparcos dPerryman et al.l 1 19971) or IVerhoelst et al.l 
20051) (Am « 4, p * 230 mas, M > 0.7M o ), we would have 
to imagine either (i) a strong dependence of the wavelength or 
(ii) an edge-on orbit with the secondary occulted by the primary. 
Both possibilities can be ruled out since (i) a differing spectral 
type would have aroused spectroscopists and (ii) an edge-on or- 
bit would have aroused people doing radial velocity measure- 
ments. Howev er, a lower mass p l anet of a few Jovian masse s, as 
propose d by llrwin et ail dl989» : lHatzes & Cochran! d 1993b and 
Brown (2007) is still a possibility. Our measurement gives an 
upper limit on its relative magnitude in the H band (Am > 7.75). 



4.3. Asymmetric brightness distribution of the stellar surface 



Radial velocities dMerlind 1 19991) a s well as photometry 
dRetter et al J 120 03: Tarr ant et al J 12007) indicate variations of a 
few days period. Photometric oscillations are especially notable, 
with amplitude variation of up to a percent, well above w hat is 
predicted by atmospheric models (Dzie mbowski et al.l200ll) . By 
putting a 1.7 x 10~ 3 1 cr upper limit on the flux of an eventual 
hotspot, our observations show that the temporal brightness os- 
cillations do not have a spatial counterpart. It means the source 
of these variations is most likely not due to convection cells 
and/or non-radial oscillations. It is interesting to note that in- 
terferometry could be a good tool to detect non-radial pulsation 
in variable stars (J3 Cephei, ...). 



4.2. Presence of a companion? 

Arcturus is often used for high-resolution spatial and spec- 
tral calibration dTuthill et al.ll2000l: iDecin et al. 2003). Such use 
makes this star both very well known and very important to 
know. This explains why, when Hipparcos flagged this star 
as a binary, it stirred an important debate i n the c ommunity. 
The absence of other observational evidence (Griffin 1998), un- 
certainties in the Hipparcos detection (Soderhielm& Mignard 
1998) and finally non-d etection with adaptive optics observa- 
tions (Tur ner et al.ll 19991) , convinced the community they could 



5. Imaging Arcturus 

5.1. The Wisard and Mira reconstruction softwares 

Model fitting confines the image to be within the range of a 
pre-defined model. This is a perfect tool to derive parameters 
of astronomical objects whose morphology is already known. 
However, it would not reveal any unexpected phenomenon, 
hence the need for less constraining image reconstruction. 

The image is sought by minimizing a so-called cost func- 
tion which is the sum of a regularization term plus data related 
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Fig. 11. Three different images of Arcturus. On the left hand, image produced using a limb-darkening power law with parameter 
from Table [3] The two others images are obtained using two different softwa res for regularized imaging. The right hand scale is 
derived using the Hipparcos parallax of 88.83 ± 0.53 mas dVan Leeuwenll2007l) . 



terms. The data terms enforce the agreement of the model image 
with the different kind of measured data (power spectrum, phase 
closures, complex visibilities, etc). The interpolation of missing 
data is allowed by the regularization and by strict constraints 
such as the positivity (which plays the role of a floating support 
constraint) and normalization. 

To validate this imaging process, we used two differ- 
ent recon struction alg orithms: Wis ard and Mira. Wisard 
(Meimon 2005; Mugnier et al. 2008) stands for "Weak-phase 
Interferometric Sample Alternating Reconstruction Device". 
Its approach consists in finding the image and the missing 
phase data jointly. T his technique is called self-c alibration in 
radio-interferometry dCornwe ll & Wilkinsonl l 19811) and has en- 
abled reliable images to be reconstructed in situations of par- 
tial phase indetermination. The strength of Wisard is that it 
combines, within a Bayesian framework, a recently developed 
noise model approx imation suited to optical interferometry data 
(Meimon et al. 2005), and an edge-preserving regularization 
(Mu gnier et al.l 2004) to deal with the sparsity of the data typ- 
ical of optical interferometry. 

Mira dThiebaut et al.l 120031) stands for "Multi-aperture 
Image Reconstruction Algorithm". Compared to Wisard, Mira 
does not explicitly manage the missing Fourier phase informa- 
tion: all missing information is handled implicitly in the data re- 
lated term. For instance, it is possible to reconstruct an image 
give n only the Fourier modulus information (power spectrum 
data; Thiebaut 2007). A second difference - in these image re- 
constructions of Arcturus - lies in the chosen prior. Instead of an 
edge preserving regularization, we used for this reconstruction 
a quadratic regularization criterion. To that end, we computed a 
prior which is a parametric model image of a stellar surface (a 
qu adratic limb darkening l aw). This method, similar to that used 
by iMonnier et al.l d2007bl) . has the particularity of requiring a 
rough model of the observed object. It will inject more informa- 
tion into the reconstruction, which in turn can give wrong results 
if the prior model is not right. 

We shall note that both Wisard and Mira algorithms can use 
various types of priors (entropy, Tikhonov, etc). It is therefore 
possible to have both algorithms using the same prior. The phase 
management - explicit or implicit - will however be different. 



The images reconstructed by Wisard and Mira are shown in 
Fig. [TT] A third representation of Arcturus is also presented. It is 
an image reconstructed from the parameters derived by fitting a 
power law limb-darkening prescription on the data (values pre- 
sented in TableO. We tentatively call such type of image recon- 
struction "Parametric imaging". Cuts of the brightness distribu- 
tions are presented Fig.Q~2] The similarity of the reconstructions 
is quite striking considering that the 2 reconstruction methods 
(i.e. data-fidelity terms), as well as the priors, are different. 



5.2. Discussion of image reconstruction 

The left hand image of Fig. [TT] shows a featureless limb- 
darkened star. It is not a surprise since the image is strictly con- 
strained by the prescription. However, the important result is the 
good fit of the prescription on our data. When doing parametric 
imaging, the x 2 is a strong information to judge the reliability 
of the image reconstruction. In this case, a reduced x 2 of 1 .9 
for 2413 degrees of freedom is a good validation of the derived 
image. 

When dealing with regularized imaging, it is more difficult 
to judge the reliability of an image reconstruction. This is be- 
cause the quantity which is minimized is a sum of a regulariza- 
tion term and a data term (generally, the x 2 )- The minimum of 
this cost function is therefore dependent on the regularization 
term, and no process is known which could use this minimum 
to judge on the quality of the reconstruction. The x 2 is still of 
interest, but does only give a partial view on the reliability of the 
reconstructed image: a reduced^ 2 close to one is important, but 
it is not a quantity by itself which will ensure the quality of the 
image reconstruction. 

Ultimately, the quality of image reconstruction will be de- 
pendent on the choice of the regularization term. The closer to 
the object the regularization term brings us, the closer to the re- 
ality our image reconstruction will be. According to this philos- 
ophy, it is important to have a good estimation of the prior. We 
recommend for optical interferometrist to use an adjustable reg- 
ularization term. Simultaneously or sequentially, a solution we 
propose consists in: (i) fit a parametric image which best de- 
scribe the data, and (ii) find the image which best fit both the 
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Fig. 12. Intensity slices through the x and y axes of the three 
images of Arcturus presented in Fig.[TT] The solid curve corre- 
spond to the parametric image, the dashed curves to the Wisard 
reconstruction, and the dotted curves to the Mira reconstruction. 
If we suppose the limb darkening prescription to be the correct 
brightness distribution of the object, we can derive the residual 
of the reconstruction obtained by the two imaging softwares: the 
rms error is around 5% for the Wisard prior (edge-preserving) 
and 2-3% for Mira (limb-darkening prior). 



data and the parametric image. This technique was the one used 
with the Mira reconstruction software. 



6. Conclusion 

In this paper, we presented an Arcturus dataset of interest- 
ing quality. With the IOTA/IONIC interferometer, we measured 
fringe contrasts of less than a percent, with errors bars in aver- 
age below that level. Using this data, we fitted several models 
and prescriptions. The closure phases were well fitted by point- 
symmetric prescriptions. No companion at less than one AU was 
detected with an upper limit on its contrast ratio of 8 x 1CT 4 . The 
same modeling of the closure phases allowed the derivation of an 
upper limit on the heterogeneity of the photosphere: no hotspot 
with a brightness above 1.7 x 10~ 3 the total flux of the photo- 
sphere was detected. 

We adjusted Marcs atmospheric models to the data. The de- 
rived Rosseland diameter equaled 21 .05 ± 0.21 mas, most of the 
error bar being induced by non-trivial wavelength calibration. 
With a reduced x 2 of 2, it is interesting to note that atmosphere 
models of regular K giants can now be challenged by interfer- 
ometry at a very fundamental level, even though spectroscopic 
agreement is near-perfect. Interestingly, we noted (i) a slight 
inconsistency in the magnitude of the limb-darkening at short 
wavelength (A « 1.55 fan; see Fig.|9]l, and (ii) a slight chromatic 
effect present in the residual (lower panel of Fig. |7}. This last re- 
sults could hint the presence of a marginal (» 0.5%) water vapor 
emission outside the photosphere. 

Finally, we imaged the photosphere using two different re- 
construction algorithms (Wisard and Mira). Both produced re- 
alistic images, but highlight the difficulty to judge the reliabil- 
ity of regularized image reconstruction. Comparatively, we pre- 
sented an image reconstructed from an ad-hoc prescription of a 
limb darkened stellar surface. The low number of free parame- 
ters, combined with a good fit to the data, hinted to us that the 
most realistic brightness distribution is in fact the one of a simple 
limb darkened disk. 
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